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We present an entanglement purification protocol for a mixture of a pure entangled state and a 
pure product state, which are orthogonal to each other. The protocol is combination of bisection 
method and one-way hashing protocol. We give recursive formula for the rate of the protocol for 
different states, i.e. the number of maximally entangled two-qubit pairs obtained with the protocol 
per a single copy of the initial state. We also calculate numerically the rate for some states. 

PACS numbers: 



Entanglement is resource in quantum information. 
Maximally entangled states are basic ingredient of funda- 
mental quantum information protocols like e.g. quantum 
teleportation pQ, dense coding [5] or Ekert's quantum 
cryptographic protocol [3J. In these protocols the max- 
imally entangled pair is shared by two parties and used 
to perform a certain task. However, in the real world the 
parties share noisy entangled pairs. Bennett et al. have 
shown that many pairs in mixed entangled states can be 
distilled to a smaller number of pairs in nearly maximally 
entangled states |U|S]. In particular, they presented pu- 
rification protocols that can be realized by means of local 
operations and classical communication (LOCC). Let us 
suppose that the parties share n pairs of qubits, each of 
which is in the state p. If using a particular protocol the 
parties can obtain m pairs of qubits, each of which is in 
the maximal entangled state, from n pairs of qubits, each 
of which is in state p then the protocol has the rate — 
for the state p. Different protocols have different rates 
and moreover one needs different protocols for different 
states. The maximal rate for state p, i.e., the rate of the 
optimal protocol for state p, is called distillable entan- 
glement. Distillable entanglement of mixed states is usu- 
ally difficult to calculate and it is only known for bound 
entangled states, maximally correlated states and some 
other specific mixed states M M EH EE] • For bound 
entangled states it is equal to zero [B], while for maxi- 
mally correlated states it is equal to the relative entropy 
of entanglement [TTJ [T2j [T3J [14] and can be distilled by 
the one-way hashing protocol |15[I16]. However, for many 
states only upper and lower bounds on distillable entan- 
glement are known. Lower bounds are usually given by 
the rates of particular protocols. Upper bounds are given 
by entanglement measures known to be greater or equal 
to distillable entanglement, e.g. the relative entropy of 
entanglement. 

In this paper we present an entanglement purification 



protocol for a mixture of a pure entangled state and a 
pure product state which are orthogonal to each other. 
The very first protocol for these states was presented in 
[5]. Our protocol is a combination of Procrustean method 
of entanglement concentration |17j , bisection method and 
one-way hashing protocol [TS] HFJ • 

We assume that Alice and Bob share many copies of 
the state 

p = p|^+)(^+| + (l-p)|00)<00|, (1) 

where 

l^ + > = ^(|10) + |01)) (2) 

is the maximal entangled state. State p is a mixture of 
the maximal entangled state and a product state. The 
product state is orthogonal to the maximal entangled 
state. Let Alice and Bob group these states in blocks of 
n copies, where n is the power of 2. It is convenient to 
write p® n in the following way 

p ®n =p «| 1 /,+ )(-0+|®« + 

p n - x {i - P )[|V' + )<V' + l® (n " 1) |oo)(oo| + ...] + 
+ p «-2(i - p) 2 [|v>+)(v> + |® ( "- 2) |oo)(oo|® 2 + . . . ] 

■•• + (l-p) n |00)(00|®", (3) 

where " . . . " in each square bracket stands for all permu- 
tations of the first term in the square bracket. 

Let each party project her/his part of the state on a 
subspace spanned by vectors with the definite number of 
l's and definite number of 0's, i.e., Alice and Bob per- 
form von Neumann measurements given by the sets of 
projectors 

{p:= e ( 4 ) 

[x]—n,\x\—a 
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and 



(5) 



[a:]— n,|a:|— 6 



respectively, where \x\ denotes the Hamming weight of 
the string x of [a;] = n qubits and a, b S {0, 1, ...n}. If 
Alice obtains P™ as a result of her measurement and 
Bob obtains P£ as a result of his measurement, then all 
terms in the expansion of p® n except the term with a + b 
|V' + )'s and n — a — b |00)'s are annihilated, i.e., the post 
measurement state is 



p(n,a,b) = 

^b)PaPb[\^ + ) |® (Q+b) 1 00) (QQf> n - a - b 



p( 



where 



p(n,a,b) 



n 
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K a + bj \ <i 
The probability of this event is 

P(n,a,b) = p a+b (l - p) n -' ■'• 



-(a+b) 



p(n,a,b). 



(6) 
(7) 
(8) 



The factor p a+5 (l - p) n_a_6 ( + 6 ) is the 
probability that Alice and Bob share the 
state |?/' + )(V' + | <3(a+&) |00)(00| 8n - a - h or any 
permutation of it and (^+ b ^2'( a + b '> 
Tr[P^P 6 n |V) + )(V' + |® (a+b) |00)(00|®"- Q - 6 P™P^]. 

If a + 6 — n, then Alice and Bob share the maximal 
entangled state of the rank r(n, a), where 



r(n, a) = 



(9) 



i.e., they share log 2 r(n,a) maximally entangled pairs of 
qubits. If one of the equalities: a = 0, a = n, & = or 
b = n holds, then Alice and Bob share a separable state. 

In the remaining cases Alice and Bob share a mixed 
entangled state. In order to distill entanglement from 
it, they can proceed in two ways. Firstly, if Alice and 
Bob share a large number of blocks of qubits in identical 
post measurement states, then they can apply the one- 
way hashing protocol [TS] and distill entanglement at the 
rates 



I C (A>B) = S(B)-S(AB), 
if Alice classicaly communicates to Bob or 
I C (B > A) = S(A) - S(AB), 



(10) 



(11) 



if Bob classically communicates to Alice. S(A) and S(B) 
are von Neumann entropies of Alice's and Bob's subsys- 
tems, respectively and S(AB) is von Neumann entropy 
of the whole system. They are given by the following 
formulae: 



S(A) = log 2 



(12) 



because the state of Alice's subsystem is an equal mixture 
of all sequences of length n with the Hamming weight 
equal to a; 



S(B) = log 2 



(13) 



because the state of Bob's subsystem is an equal mixture 
of all sequences of length n with the Hamming weight 
equal to b; 



S(AB) = log 2 



(14) 



because the state of the whole system is an equal mixture 
of ( a " b ) pure orthogonal states. Hence, the optimal rate 
of the one-way hashing protocol is 



log 2 



I c (n,a,b) - 



log 2 



i—a—b/ 



(15) 



Secondly, Alice and Bob can divide the pairs of qubits 
into two blocks of equal length (the first block consists of 
the first n/2 pairs and the second block consists of the last 
n/2 pairs) and repeat the measurements on each block 
separately. If Alice obtains P™, 
her measurement and Bob obtains Py 



P™/ 2 as the result of 



,n/2 



P fc ™/ 2 as the 



result of his measurement, then the post measurement 
state is 



p(n/2,a',b') ® p(n/2,a" ,b"), 



(16) 



where a' + a' 1 
the identity 



a and b' + b" — b. In derivation we used 



(p;/ 2 ® p n J 2 )p: = (p 7 :/ 2 ® p;/ 2 )<W',*- (17) 



The probability of this event is 



u jj „ m v\ P(n/2,a',b')p(n/2,a",b") 

p(a,b;a ,b \n,a,b) = .(18) 

p(n,a,b) 



If a' + b' = n/2 (a" + b" = n/2), then the first (sec- 
ond) block of pairs is in the maximal entangled state 
of the rank r(n/2,af) (r(n/2, a"))- If a' + b' ^ n/2 
(a" + b" 7^ n/2), then Alice and Bob can choose if they 
want to apply the one-way hashing protocol or to di- 
vide the pairs in the first (second) block into two blocks 
and repeat the measurement on each block separately. 
For different choices Alice and Bob obtain different rates. 
The rates achievable with the optimal choices are given 
by the following recursive formula: 
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R(n, a, b) = max{/ c (n, a, b), 
Eab^ Eht,> mixi P(a', b'; a-a',b- b'\n, a, b)(R(n/2, a', b>) + R(n/2, a-a',b- &'))} (19) 



where the summation limits are 

a min = max{0,a- §} 
a max = min{a, f } 
b' min = max{0, a + b — a' — ^} 

b' max = min{6, f - a'} (20) 

If we use only the bisection method (without the one- 
way hashing protocol), then we obtain the following ex- 
pression for the rate of such a protocol for state p: 

m 

R{ P ) = Y,p {2l) W)-y^ l -% (2i) 
i=i 

where n — 2 m and 

In Fig. [l]we present the rates of the protocol based on 
the bisection method and the one-way hashing protocol 
for different initial states p of Eq. [T] which depend on 
the parameter p. The first measurement was performed 
on a block of 64 pairs of qubits. For comparison, we 
also present the rates of the protocol from [5] and the 
one-way hashing protocol as well as an upper bound for 
distillablc entanglement given by the relative entropy of 
entanglement. One can see that our protocol performs 
better than the other two protocols. 



R 




FIG. 1: Relation between coherent information (thin solid 
line), relative entropy of entanglement (thick solid line), the 
rate of the protocol of Bennett et al. (dashed line) and our 
protocol for n — 64 (dots) for different values of the parameter 
P- 



measurement was performed. For comparison, we present 
also the rates of the protocol based only on the bisection 
method (without the one-way hashing protocol). One 
can see that the one-way hashing protocol causes an in- 
crease in the rate. 

TABLE I: Results for p = |. n - size of the block of pairs of 
qubits on which the first measurement is performed. 7? - the 
rate of the protocol based on bisection method and one-way 
hashing protocol. R' - the rate of the protocol based only on 
bisection method. 



n 


R 


R' 


2 


0.111111 


0.111111 


4 


0.158981 


0.158981 


8 


0.173419 


0.16638 


16 


0.175076 


0.166574 


32 


0.175129 


0.166575 


64 


0.175129 


0.166575 



Let us also point out that our protocol applies to a 
mixture of an arbitrary pure entangled state and a pure 
product state orthogonal to it, i.e., 

p = p|V)(^| + (l-p)|00)<00| (23) 

where 

|V) = a|10)+/3|01) (24) 

One can see it by noting that the first measurement 
projects n copies of the initial state on a state given in 
Eq. [6] with probability 

P'(n,a,b) =p a + b (l - P ) n - a - b ( a l b ) x 

x( a :>H/?r (25) 

Because the states given in Eq. [23] can be obtained by 
sending half of the state 

|V/) = a'|10) +/ 3'|01)) (26) 

through the amplitude damping channel Nad, given by 
Kraus operators 



In Table [T] we present the rates of the protocol for dif- 
ferent sizes of blocks of pairs of qubits on which the first 



S = |0)(0| + v ^f|l)(l| 

£i= 7 |0)(l| (27) 
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we can obtain the lower bound on quantum capacity of 
the amplitude damping channel assisted by two-way clas- 
sical communication Q2- In such a case, the best lower 
bound on Q2 is given by the following expression 

R mx = mBxR(N AD (il/)). (28) 

i/i' 

Our protocol requires collective measurements on large 
blocks of pairs of qubits, while the protocol presented in 
[5] requires only measurements on two pairs of qubits. 
We will show how one can improve the two-copy protocol. 
Let Alice and Bob perform von Neumann measurements 
given by projectors 

P = |Q0)(00| + |11)(11| 

P 1 = |01>(01| + |10)(10| (29) 

on the state p® 2 If both obtain 1 as the results of their 
measurements, then the post measurement state is equiv- 
alent to the maximal entangled state \ip + ). The proba- 

2 

bilty of this event is equal to If both obtain as the 
results of their measurements, then the post measure- 
ment state is equivalent to 

p'=p>+)(V' + | + (l-p')|00)<00| (30) 

2 

where p' — p 2 + 2(i~p) 2 • ^he probability of this event is 

equal to \ + (1 — p) 2 . Two such states can be used 
in another measurement and hence there is a chance to 
obtain from them the maximal entangled state. In par- 
ticular, if p = 2/3 then p' = 2/3 and the rate is equal to 
2/15 for improved protocol instead of 1 /9 for the original 
protocol. 

In conclusion, we presented the entanglement purifica- 
tion protocol for a mixture of a pure entangled state and 
a pure product state, orthogonal to each other. To our 
knowledge uur protocol performs better than any other 
previously known protocol. We also discussed how one 
can obtain the lower bound on quantum capacity of the 
amplitude damping channel assisted by two-way classical 
communication. 

Fruitful discussion with M. Horodecki is acknowledged. 
One of the authors (A.G.) was partially supported by 
Ministry of Science and Higher Education Grant No. N 
N206 2701 33 and by the European Commission through 
the Integrated Project FET/QIPC SCALA. 



[4] C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, 
J. A. Smolin, and W. K. Wootters, Phys. Rev. Lett. 76, 
722 (1996). 

[5] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, and W. K. 

Wootters, Phys. Rev. A 54, 3824 (1996). 
[6] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 

Rev. Lett. 80, 5239 (1998). 
[7] E. M. Rains, arXiv:quant-ph/9707002 (1997). 
[8] J. Eisert, T. Felbinger, P. Papadopoulos, M. B. Plenio, 

and M. Wilkens, Phys. Rev. Lett. 84, 1611 (2000). 
[9] Y.-X. Chen and D. Yang, Phys. Rev. A 66, 014303 

(2002). 

[10] S. Hamieh and H. Zaraket, J. Phys. A: Math. Gen. 36, 
L387 (2003). 

[11] T. Hiroshima and M. Hayashi, Phys. Rev. A 70, 030302 
(2004). 

[12] V. Vedral, M. B. Plenio, M. A. Rippin, and P. L. Knight, 

Phys. Rev. Lett. 78, 2275 (1997). 
[13] V. Vedral and M. B. Plenio, Phys. Rev. A 57, 1619 

(1998). 

[14] E. M. Rains, IEEE Trans. Inf. Theory 47, 2921 (2001). 
[15] I. Devetak and A. Winter, Phys. Rev. Lett. 93, 080501 
(2004). 

[16] I. Devetak and A. Winter, Proc. R. Soc. Lond. A 461, 
207 (2005). 

[17] C. H. Bennett, H. J. Bernstein, S. Popescu, and B. Schu- 
macher, Phys. Rev. A 53, 2046 (1996). 



[1] C. H. Bennett, G. Brassard, C. Crepeau, R. Jozsa, 
A. Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 
(1993). 

[2] C. H. Bennett and S. J. Wiesner, Phys. Rev. Lett. 69, 
2881 (1992). 

[3] A. K. Ekert, Phys. Rev. Lett. 67, 661 (1991). 



